Chapter 2 

Motion of Charged Particles in Fields 



Plasmas are complicated because motions of electrons and ions are determined by the electric 
and magnetic fields but also change the fields by the currents they carry. 

For now we shall ignore the second part of the problem and assume that Fields are Prescribed. 

Even so, calculating the motion of a charged particle can be quite hard. 

Equation of motion: 

m% = q ( E + v A B ) (2.1) 

^_Jt, charge E-field velocity B-field 

Rate of change of momentum ' Lorentz Force ' 

Have to solve this differential equation, to get position r and velocity (v= r) given E(r, t), B(r, t). 
Approach: Start simple, gradually generalize. 

2.1 Uniform B field, E = 0. 

mv = qw A B (2.2) 

2.1.1 Qualitatively 

in the plane perpendicular to B: Accel, is perp to v so particle moves in a circle whose 
radius r L is such as to satisfy 

v 2 

mr L tt 2 = m— = \q\v±B (2.3) 
'7. 

Q is the angular (velocity) frequency 

1st equality shows f2 2 = v\jr\ (r L = v±/Q) 
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Figure 2.1: Circular orbit in uniform magnetic field. 



Hence second gives m^fr = \q\v±B 



m 

Particle moves in a circular orbit with 



i.e. Q = . (2.4) 



angular velocity Q = the "Cyclotron Frequency" (2-5) 

m 

v± 

and radius r { = — the "Larmor Radius. (2.6) 

2.1.2 By Vector Algebra 

• Particle Energy is constant, proof: take v. Eq. of motion then 

mv.v = — (^-mv 2 ^j = qv.(v A B) = 0. (2.7) 

• Parallel and Perpendicular motions separate. v\\ = constant because accel (oc v A B) 
is perpendicular to B. 

Perpendicular Dynamics: 

Take B in z direction and write components 

mv x = qv y B , mv y = —qv x B (2.8) 

Hence 



qB . (qB 
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v x = — i)y = - — v x — —Q v x (2.9) 
m \m ' 



Solution: v x = v± cos Qt (choose zero of time) 
Substitute back: 



ttl q 

Vy = -j^Vx = v±sinQt (2-10) 
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Integrate: 



x = 



v±_ . „ q v± 

s + — smSH , y = Ho ~\~ — r cos j it 



(2.11) 




q positive (ion) 
Left-handed 
Anticlockwise 

q negative (e) 
Right-handed 
Clockwise 



Figure 2.2: Gyro center (xo,yo) and orbit 

This is the equation of a circle with center r = (x , y ) and radius r L = v±/Q: Gyro Radius. 
[Angle is 6 = fit] 

Direction of rotation is as indicated opposite for opposite sign of charge: 

Ions rotate anticlockwise. Electrons clockwise about the magnetic field. 

The current carried by the plasma always is in such a direction as to reduce the magnetic 
field. 

This is the property of a magnetic material which is "Diagmagnetic" . 
When v\\ is non-zero the total motion is along a helix. 



2.2 Uniform B and non-zero E 

mv = q(E + v A B) 
Parallel motion: Before, when E = this was v\\ = const. Now it is clearly 

v... - qE W 



m 



(2.12) 



(2.13) 



Constant acceleration along the field. 



Perpendicular Motion 
Qualitatively: 

Speed of positive particle is greater at top than bottom so radius of curvature is greater. 
Result is that guiding center moves perpendicular to both E and B. It 'drifts' across the 
field. 

Algebraically: It is clear that if we can find a constant velocity that satisfies 



E + v d A B = 



(2.14) 
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E 



B 




Figure 2.3: E A B drift orbit 



then the sum of this drift velocity plus the velocity 

d , 



which we calculated for the E = gyration will satisfy the equation of motion. 
Take AB the above equation: 

= E A B + (v<j AB)AB = EAB + (u d .B)B - B 2 v d 



so that 

does satisfy it. 

Hence the full solution is 



E AB 

B 2 



where 



and 



V = V|| + v d + V L 

parallel cross-field drift Gyration 



1 E \\ 



rn 



(2.15) 



(2.16) 
(2.17) 



(2.18) 



(2.19) 



\d (eq |2.17 ) is the "E x B drift" of the gyrocenter. 
Comments on E x B drift: 

1. It is independent of the properties of the drifting particle (q, m, v, whatever). 

2. Hence it is in the same direction for electrons and ions. 

3. Underlying physics for this is that in the frame moving at the E x B drift E = 0. We 
have 'transformed away' the electric field. 

4. Formula given above is exact except for the fact that relativistic effects have been 
ignored. They would be important if Vd ~ c. 
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2.2.1 Drift due to Gravity or other Forces 

Suppose particle is subject to some other force, such as gravity. Write it F so that 

mv = F + gvAB = g(VvAB) (2.20) 

This is just like the Electric field case except with F/q replacing E. 
The drift is therefore 

1 F A B 

v^-p- (2.21) 

In this case, if force on electrons and ions is same, they drift in opposite directions. 

This general formula can be used to get the drift velocity in some other cases of interest (see 
later) . 

2.3 Non-Uniform B Field 

If B-lines are straight but the magnitude of B varies in space we get orbits that look quali- 
tatively similar to the E 1 B case: 



VB 




Figure 2.4: V-B drift orbit 

Curvature of orbit is greater where B is greater causing loop to be small on that side. Result 
is a drift perpendicular to both B and VB. Notice, though, that electrons and ions go in 
opposite directions (unlike E A B). 

Algebra 

We try to find a decomposition of the velocity as before into v = + v L where v rf is 
constant. 

We shall find that this can be done only approximately. Also we must have a simple expres- 
sion for B. This we get by assuming that the Larmor radius is much smaller than the scale 
length of B variation i.e., 

r L «B/\VB\ (2.22) 
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in which case we can express the field approximately as the first two terms in a Taylor 
expression: 

B ~ B + (r.V)B (2.23) 
Then substituting the decomposed velocity we get: 

= mv t = ?(vAB) = g[v L A B + v d A B + (v L + v rf ) A (r.V)B] (2.24) 
or = v d AB + v L A (r.V)B + v d A (r.V)B (2.25) 

Now we shall find that Vd/vi is also small, like r\VB\/B. Therefore the last term here is 
second order but the first two are first order. So we drop the last term. 

Now the awkward part is that Vl and r L are periodic. Substitute for r = r + r L so 

= w d A B + v L A (r L .V)B + v L A (r .V)B (2.26) 
We now average over a cyclotron period. The last term is oc e~ lQ,t so it averages to zero: 

= v rf A B + (v L A (r L .V)B) . (2.27) 

To perform the average use 



r L = (xL,y L ) = ^ (^mfit^cosfitj (2.28) 
v L = (x L ,y L ) = v± (cos fit, sin fit j (2.29) 



n 


^sinfit, 


\q\ 


v± 


^cos fit, 


\q\ 


v y y 


dy B 





So [vl A (r.V)B] s = v y y—B (2.30) 



[v L A(r.V)B] y = -v x y^B (2.31) 



(Taking VB to be in the y-direction). 
Then 



v 2 



So 



Substitute in: 



( Vy y) = - (cos fit sin fit) -± = (2.32) 

2 j 2 

(v x y) = (cos fit cos fit) ^ = 2^y^j ( 2 - 33 ) 



(v L A(r.V)B) = -j|^VS (2.34) 



= v d AB-±£vB (2.35) 
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and solve as before to get 



^VB AB „ „_> 



Vd ~ W -\q\m^^' (2 ' 36) 



or equivalently 

_ lm^ BAVB 

~~ g~2£? 52 (2 ' 37) 

This is called the 'Grad B drift'. 



2.4 Curvature Drift 

When the B-field lines are curved and the particle has a velocity v\\ along the field, another 
drift occurs. 




Center of Curvature 

Figure 2.5: Curvature and Centrifugal Force 

Take \B\ constant; radius of curvature R e . 

To 1st order the particle just spirals along the field. 

In the frame of the guiding center a force appears because the plasma is rotating about the 
center of curvature. 

This centrifugal force is F c f 

II 

F c f = 771 -zj- pointing outward (2.38) 



as a vector 



F cf = mv lj£ ( 2 - 39 ) 



[There is also a coriolis force 2m(w A v) but this averages to zero over a gyroperiod.] 
Use the previous formula for a force 

_ 1 F cf A B mvj R e A B 
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This is the "Curvature Drift". 

It is often convenient to have this expressed in terms of the field gradients. So we relate R c 
to V.B etc. as follows: 




Figure 2.6: Differential expression of curvature 



(Carets denote unit vectors) 
From the diagram 

and 

So 

But (by definition) 



db = b 2 - bi = -R c tt (2.41) 
di = oc R c (2.42) 



dX) R c R c 

~dl = ~R~ C = ~B? C 



dl 

So the curvature drift can be written 

,2 T, t-> ™„,2 



2.4.1 Vacuum Fields 



Relation between VB & R r drifts 



(2.43) 



' /b (B.V)b (2.44) 



_ m^R c B mt, B A (b,V)b 

Vd "^if A ^-^ ^ (2 ' 45) 



The curvature and V-B are related because of Maxwell's equations, their relation depends 
on the current density j. A particular case of interest is j = 0: vacuum fields. 



B w 



r{= R') 

Figure 2.7: Local polar coordinates in a vacuum field 

V A B = (static case) (2.46) 
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Consider the z-component 



1 d 

= (V A B) g = - — {rB e ) (B r = by choice). (2.47) 
r or 

° D " + * (2.48) 



dr r 

or, in other words, 



(VB) r = ~ (2.49) 



[Note also = (V A B) e = dB e /dz : (VB) Z = 0] 
and hence (V-B) pcrp = — B ^R C /R 2 C . 
Thus the grad B drift can be written: 



_ mv\ B A VB _ mVl R c A B 
VvB "^g W~~ 2q R 2 C B 2 ^ 



and the total drift across a vacuum field becomes 



1 ( 2 1 2 A R c A B 



Vr + v V b = - (m^, + -mv ± j R2B2 . (2.51) 

Notice the following: 

1. R c & Vi? drifts are in the same direction. 

2. They are in opposite directions for opposite charges. 

3. They are proportional to particle energies 

4. Curvature <-> Parallel Energy (x 2) 

V-B <-> Perpendicular Energy 

5. As a result one can very quickly calculate the average drift over a thermal distribution 
of particles because 

{\mv\) = T - (2.52) 
(-mvD = T 2 degrees of freedom (2.53) 



Therefore 



"2 



, 2TR C AB / 2TB A (b.Vjb 

(v fi + v TO ^-— U- ^ 1 (2.54) 
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2.5 Interlude: Toroidal Confinement of Single Parti- 
cles 



Since particles can move freely along a magnetic field even if not across it, we cannot ob- 
viously confine the particles in a straight magnetic field. Obvious idea: bend the field lines 
into circles so that they have no ends. 




Figure 2.8: Toroidal field geometry 



Problem 

Curvature & VB drifts 



If 2 1 2 \ RAB 

1 ( 2 1 2 \ 1 



(2.55) 
(2.56) 



Ions drift up. Electrons down. There is no confinement. When there is finite density things 




EAB 



Figure 2.9: Charge separation due to vertical drift 
are even worse because charge separation occurs — > E — > E A B — > Outward Motion. 



2.5.1 How to solve this problem? 

Consider a beam of electrons ^ fj_ = 0. Drift is 

.2 T 



Vd = 



q B T R 
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(2.57) 



What B z is required to cancel this? 

Adding B z gives a compensating vertical velocity 

B 7 



We want total 



v = uti — — for B z << B T 

Bj< 



B z mv\ q 

v z = = v\\ 1 

Bt q Bj'R 



(2.58) 



(2.59) 



So B z = —mv\\/Rq is the right amount of field. 
Note that this is such as to make 



r L {B z 



\mvu | 



R 



(2.60) 



But B z required depends on v\\ and q so we can't compensate for all particles simultaneously. 
Vertical field along cannot do it. 

2.5.2 The Solution: Rotational Transform 




Figure 2.10: Tokamak field lines with rotational transform 
Toroidal Coordinate system (r, 9, (p) (minor radius, poloidal angle, toroidal angle), see figure 



Suppose we have a poloidal field Bg 



Field Lines become helical and wind around the torus: figure 2.10 
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In the poloidal cross-section the field describes a circle as it goes round in 4>. 
Equation of motion of a particle exactly following the field is: 

Bn 



d9 B d Bq Bm 



B 



and 



r = constant. 

Now add on to this motion the cross field drift in the z direction. 




Figure 2.11: Components of velocity 



dO 

di 
dr 

dt 



Be 



v\\ + Vd cos 9 



B 

Vd sin 9 



Take ratio, to eliminate time: 



1 dr Ud sin 9 

r d9 %Lv\\ + v d cos 9 

Take Be, B,v\\,Vd to be constants, then we can integrate this orbit equation: 

[lnr] = [-ln|^i + t; d cosfl|] . 



Take r = r when cos^ = (9 = |) then 



r = r / 



i , Bvd a 
1 + cos 9 

Mil 



If J^j- << 1 this is approximately 



r = r — A cos 9 



r 



where A = J^- 
This is approximately a circular orbit shifted by a distance A: 
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Figure 2.12: Shifted, approximately circular orbit 



Substitute for v d 



B 1 mt) + \mv\) 1 

r o^-- — 2.69 

B e q v\\ B^R V ; 

1 mv\ + \mvlr p 



qB e v\\ R 

jf n A mv \\ r ° r ° 71 ^ 

UV± = ° A = WeR =rL6 R> (2 - 71) 

where r^g is the Larmor Radius in a field Bg x r/R. 

Provided A is small, particles will be confined. Obviously the important thing is the poloidal 
rotation of the field lines: Rotational Transform. 

Rotational Transform 

poloidal angle , n „. 

rotational transform = — (2.72) 

1 toroidal rotation 

(transform/27r =) t= P°^°^ a | an g| e ^ (2.73) 

toroidal angle 

(Originally, t was used to denote the transform. Since about 1990 it has been used to denote 
the transform divided by 2n which is the inverse of the safety factor.) 

'Safety Factor' 

, , 1 toroidal angle , _ , 

V a = - = — r^n — V • 2 - 74 

i poloidal angle 

Actually the value of these ratios may vary as one moves around the magnetic field. Definition 
strictly requires one should take the limit of a large no. of rotations. 

q s is a topological number: number of rotations the long way per rotation the short way. 

Cylindrical approx.: 
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In terms of safety factor the orbit shift can be written 



I A I r B ^ T 

|A| = r L6 — = r L4> —— = r L q s 



(2- 



(assuming B^ » B e ) 



2.6 The Mirror Effect of Parallel Field Gradients: E 

0, VB II B 




Figure 2.13: Basis of parallel mirror force 



In the above situation there is a net force along B. 
Force is 



< F|| > = — |gvAB| sina;= — \q\v±B sin a 
-B r 



sin a = 



B 



Calculate B r as function of B z from V.B = 0. 



Hence 



1 d d 
V.B = - — (rB r ) + —B z = 0. 
r or oz 



rB r = — I r—^-dr 



dz 



Suppose 7"l is small enough that ~ const 

[rB r ]? ~ [ 
Jo 



rL , dB z 
rdr — — 
oz 



1 2 dB z 

2 T ' L dz 



So 



r> / \ 1 dB z 



B r 



r L 1 dB z 



since = — — = +- 

B 2 2 dz 



(2- 
(2- 

(2- 
(2- 

(2- 

(2- 
(2- 
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Hence 

^ p ^_ \ n \V±r L dB g _ \mv\ dB z 

As particle enters increasing field region it experiences a net parallel retarding force. 
Define Magnetic Moment 

/I = \mv\jB . (2.85) 
Note this is consistent with loop current definition 

AT 2 l?K WlV±_ . , 

ii — AI — nr T . = 2.86 

Force is F\\ = //.VyB 

This is force on a 'magnetic dipole' of moment fj,. 

F|| = //.V||B (2.87) 

Our always points along B but in opposite direction. 

2.6.1 Force on an Elementary Magnetic Moment Circuit 

Consider a plane rectangular circuit carrying current I having elementary area alxaly = dA. 
Regard this as a vector pointing in the z direction dA. The force on this circuit in a field 
B(r) is F such that 

OB 

F x = Idy[B z (x + dx) - B z (x)\ = Idydx—^- (2.88) 

OB 

F y = -Idx[B z (y + dy)-B z (y)]=Idydx-^- (2.89) 
F z = -Idx[B y (y + dy)-B y (y)]-Idy[B x (x + dx)-B x (x)} (2.90) 
= —Idxdy 



dB x dB v 



dx dy 



Idydx—^- (2.91) 

az 



(Using V.B = 0). 

Hence, summarizing: F = IdydxV B z . Now define \i = IdA = Idydxz and take it constant. 
Then clearly the force can be written 

F = V(B.fi) [Strictly = (VB)./x] (2.92) 
H is the (vector) magnetic moment of the circuit. 

The shape of the circuit does not matter since any circuit can be considered to be composed 
of the sum of many rectangular circuits. So in general 

/i = IdA (2.93) 
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and force is 

F = V(B./i) (fi constant), (2.94) 

We shall show in a moment that is constant for a circulating particle, regard as an 
elementary circuit. Also, fi for a particle always points in the -B direction. [Note that this 
means that the effect of particles on the field is to decrease it.] Hence the force may be 
written 

F = -fi\7B (2.95) 



This gives us both: 

• Magnetic Mirror Force: 

and 

• Grad B Drift: 



F\\ = -nV\\B (2.96) 



1 FAB /iBAVfl , on _, 



2.6.2 fi is a constant of the motion 

'Adiabatic Invariant' 



Proof from Fu 

Parallel equation of motion 



dv\\ dB 



So 



or 



dvu d ,1 2 \ dB dB 

mv\\—r- = -r(-mvu) = —uv z —— = —u—r- (2.99) 
11 dt dt K 2 11 ; p dz p dt K ' 

d . 1 on dB 



Conservation of Total KE 



Combine 



Y{\ mv l + \ mv l) = (2.101) 

(XL Zi Zi 

^-mv\ + siB) = a (2.102) 



= ^ = As required (2.104) 
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Angular Momentum 



of particle about the guiding center is 

mv± 2m \mv\ 

r L mv± = T - rF rmv± = — r „ (2.105) 
\q\B \q\ B 

= ^/i . (2.106) 

m 

Conservation of magnetic moment is basically conservation of angular momentum about the 
guiding center. 

Proof direct from Angular Momentum 

Consider angular momentum about G.C. Because 6 is ignorable (locally) Canonical angular 
momentum is conserved. 

p = [r A (mv + qA)] z conserved. (2.107) 

Here A is the vector potential such that B = V A A 
the definition of the vector potential means that 

ft = 1 (2.08) 
r or 



T 2 

r L A e (r L ) = [ L r.B z dr= r ^B z = ^ (2.109) 
Jo 2 \q\ 



Hence 



— q mfi 

p = -j— 7-rLV±m + q— — (2.110) 



in \q 
\q 



-T^m/x. (2.111) 



So p = const = constant. 

Conservation of /x is basically conservation of angular momentum of particle about G.C. 
2.6.3 Mirror Trapping 

F|| may be enough to reflect particles back. But may not! 
Let's calculate whether it will: 
Suppose reflection occurs. 
At reflection point v\\ r = 0. 
Energy conservation 

2™( u lo + 4>) = 2 mu lr (2-112) 
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Reflection 




B-lines 



Stronger 



pi conservation 



Hence 



Figure 2.14: Magnetic Mirror 



k mv ±o 



\mv\ r 



Bo 



B r 



2 , 2 
V ±0 + V \\0 = 


B r 2 
Bo 


Bo 




B r 





2.6.4 Pitch Angle 9 



tan 9 

Bo 

B, r 



v± 

V\\ 



v lo + v \\o 



sin 2 6*o 



So, given a pitch angle 9q, reflection takes place where B /B r = sin 2 9q. 
If 6q is too small no reflection can occur. 
Critical angle 9 C is obviously 

9 C = sin' 1 (B /B 1 )^ 

Loss Cone is all 9 < 9 C . 

Importance of Mirror Ratio: R m = Bi/B . 



2.6.5 Other Features of Mirror Motions 

Flux enclosed by gyro orbit is constant. 



2 2 

■Km v , 



$ = uriB = 

q 2 B 2 
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B 



Figure 2.15: Critical angle 9 C divides velocity space into a loss-cone and a region of mirror- 
trapping 



9 

Note that if B changes 'suddenly' ji might not be conserved 



MS 



Flux Tube 




(2.120) 



2-nm \mv\ 
2nm 

— —[/, = constant. (2.121) 



Figure 2.16: Flux tube described by orbit 

Basic requirement 

r L « B/\VB\ (2.122) 

Slow variation of B (relative to r L ). 

2.7 Time Varying B Field (E inductive) 

Particle can gain energy from the inductive E field 

<9B 

V A E = -— (2.123) 
or j E.dl = -^B.ds = -^ (2.124) 

Hence work done on particle in 1 revolution is 

/( • d<& 
\q\E.d£ = +\q\ J B.ds = +\<l\^ = \q\Bmj, (2.125) 
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Figure 2.17: Particle orbits round B so as to perform a line integral of the Electric field 
(d£ and vj_g are in opposition directions). 

= 1^-1 = ^^ (2-126) 



Hence 



but also 



Hence 



d (\ 9 \ M /l ,\ cfc 

sM-krHa™ 4 )""* (2 - 128) 



|(^) = |(MB). (2.29) 
| = 0. (2.130) 



Notice that since $ = ^pr/i, this is just another way of saying that the flux through the gyro 
orbit is conserved. 

Notice also energy increase. Method of 'heating'. Adiabatic Compression. 



2.8 Time Varying E-field (E, B uniform) 

Recall the E A B drift: 



EAB 

when E varies so does veab- Thus the guiding centre experiences an acceleration 



VEAB = (2.131) 



d /EAB\ 



In the frame of the guiding centre which is accelerating, a force is felt. 

d /EAB 

dt I - B 



F a = —m— ( — — — J (Pushed back into seat! — ve.) (2.133) 
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This force produces another drift 



1 F a A B m d /EAB\ 

= - — ^f(E.B)B- J B 2 E) 

qB dt^ J > 

m ■ 



qB 2 



This is called the 'polarization drift 7 . 



(2.134) 

(2.135) 
(2.136) 



Vfl = VE/\B + Vp = 



E AB 



m 



B 2 qB 2 
E AB 1 



B 2 



QB 



+ — E 
E 



(2.137) 
(2.138) 




Initiar^osition at rest (when E = 0). Final average (gyro) locus. 



Figure 2.18: Suddenly turning on an electric field causes a shift of the gyrocenter in the 
direction of force. This is the polarization drift. 

Start-up effect: When we 'switch on' an electric field the average position (gyro center) of 
an initially stationary particle shifts over by ~ \ the orbit size. The polarization drift is this 
polarization effect on the medium. 

Total shift due to v p is 



Ar / Vpdt 



m 



qB 2 



E i dt 



m 



qB 2 



[AE^ 



(2.139) 



2.8.1 Direct Derivation of ^ effect: 'Polarization Drift 



dt 



Consider an oscillatory field E = Ee twt 

dv 



m 



dt 



(± r B) 

q (E + v A B) 
q (Ee" 



+ v A 



Try for a solution in the form 



v = v D e-^' + v L 



(2.140) 
(2.141) 

(2.142) 
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where, as usual, v L satisfies mw L = q\ L A B 
Then 



-iuit 



(1) m(-icuv D = q (E + v D A B) xt 
Solve for v D : Take AB this equation: 

(2) - miuj (v D A B) = q (E A B + (B 2 .v|£>) B - B 2 v D 
add miuj x (f) to q x (2) to eliminate v D A B. 

m 2 uj\ D + q 2 (E A B - £ 2 vd) = miuoqE 

E AB 



or : 



i.e. 



Vfl 



f - 



2 2 

q 2 B 2 

n 2 



miu 
iuiq 



E + 



fi 2 
E AB 



(2.143) 
(2.144) 

(2.145) 
(2.146) 

(2.147) 



Since — iu 



dt 



this is the same formula as we had before: the sum of polarization and 



E A B drifts except for the [1 — uj 2 VL 2 } term. 
This term comes from the change in with time (accel). 
Thus our earlier expression was only approximate. A good approx if uj « Q. 



2.9 Non Uniform E (Finite Larmor Radius) 



m 



dt 



q (E(r) + v A B) 



Seek the usual soltuion v = v D + v 9 . 
Then average out over a gyro orbit 

dv 



Hence drift is obviously 



(m^) = = (q (E(r) + v A B)) 
= g[(E(r))+v D AB] 

(E(r))AB 
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So we just need to find the average E field experienced. 
Expand E as a Taylor series about the G.C. 

/ x 2 d 2 v 2 d 2 \ 
E(r) = E + (r.V) E + ( — + | — J E + cross terms + 



(2.148) 



(2.149) 
(2.150) 

(2.151) 



(2.152) 
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(E.g. cross terms are xy^ y 'E). 
Average over a gyro orbit: r = r L (cos9, sin6*, 0). 
Average of cross terms = 0. 
Then 



(E(r))=E + ((r L >.V)E + ^V 2 E. ( 



linear term (n) = 0. So 



2 

(E(r)) ~ E + T -±V 2 E ( 



Hence E A B with 1st finite-Larmor-radius correction is 

/ r| „\EAB 
^AB=(^l + fV 2 J^^. ( 

[Note: Grad B drift is a finite Larmor effect already.] 
Second and Third Adiabatic Invariants 

There are additional approximately conserved quantities like \i in some geometries. 



2.10 Summary of Drifts 



E A B 

v £ — — Electric Field 
B 2 

1FAB n IP 

v F = — — General Force 

q B 2 

(j.2 \ E A B 
1 + — V 2 Nonuniform E 
4 ) B 2 

mv\ B A VB _ _ 

vvb = -z ™ — GradB 

2q B 6 

mvf R c A B 
vr = „ 9 „ 9 Curvature 

1 / 2 1 2 \R C AB „ , , 

Vfl + vvB = - rra;,, + -mu 1 n9 Vacuum Fields, 

g V 11 2 / /i 2 !* 2 

o E_\_ 

v p = — Polarization 



\q\ \n\B 



Mirror Motion 



Force is F = —jiVB. 



\i = is constant ( 

2B 
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